New Experimental Limit on the Electric Dipole Moment of the Electron in a 

Paramagnetic Insulator 



O 

(N 

Sh . 
Ph. 

< 

(N 
(N 



X 
=5 



On 

o 



Y. J. Kim 1 !] C.-Y. Liu 1 !] S. K. Lamoreaux 2 , G. Visser 1 , B. Kunkler 1 , A. V. Matlashov 3 , and T. G. Reddy 1 
1 CEEM, Physics Department, Indiana University, Bloomington, IN 47408, USA 
2 Physics Department, Yale University, New Haven, CT 06520, USA and 
3 P-21, Los Alamos National Laboratory, New Mexico, NM 87544, USA 
(Dated: January 20, 2013) 

We report results of an experimental search for the intrinsic Electric Dipole Moment (EDM) 
of the electron using a solid-state technique. The experiment employs a paramagnetic, insulating 
gadolinium gallium garnet (GGG) that has a large magnetic response at low temperatures. The 
presence of the eEDM would lead to a small but non-zero magnetization as the GGG sample is 
subject to a strong electric field. We search for the resulting Stark-induced magnetization with a 
sensitive magnetometer. Recent progress on the suppression of several sources of background allows 
the experiment to run free of spurious signals at the level of the statistical uncertainties. We report 
our first limit on the eEDM of (-5.57 ± 7.98 ± 0.12) x 10~ 25 e-cm with 5 days of data averaging. 



I. INTRODUCTION 

The search for the electric dipole moment (EDM) of el- 
ementary particles is motivated to test the discrete sym- 
metries assumed in the Standard Model (SM) of par- 
ticle physics. An account of the different transforma- 
tion properties of the EDM (a polar vector) and the 
spin (a pseudo- vector), the fundamental physical laws 
governing particles must violate both time-reversal (T) 
and parity-inversion (P) symmetries for a fermion to ac- 
quire an EDM [l[ . While the phenomena of P violation 
is firmly established in numerous experiments, T viola- 
tion has only been observed directly in the neutral kaon 
system 0], with more searches in the B system under- 
way [|| . Measurements of EDMs of elementary particles 
use different experimental techniques, often on low en- 
ergy systems at Q ~ 0, to probe the physics of T vio- 
lation (with no flavor-changing) at energy scales higher 
than tens of TcV, and could provide information comple- 
mentary to that from high-energy collider experiments 
on the nature of symmetry breaking. 

The physics of T violation is often linked, via the 
CPT theorem, to the violation of the combined Charge 
conjugate (C) and P symmetry. The only confirmed 
source of CP violation in the SM is the complex phase 
of the CKM matrix (that describes the quark mixing in 
charged-current weak interactions). With it, the elec- 
tron EDM (eEDM) manifests through high-order loop 
couplings that involve flavor-changing quark interactions 
with the exchange of weak bosons. The resulting 
size of eEDM predicted within the framework of SM is 
no larger than 10 _38 e-cm, which is well beyond the reach 
of current experimental techniques. The current experi- 
mental upper bound is established using thallium atomic 
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the SM often lead to a sizable eEDM. Free from the SM 
backgrounds, the measurement of EDM presents a power- 
ful tool for global tests on many theoretical extensions to 
the SM. In particular, some variants of the popular super- 
symmetric model can generate the EDM of elementary 
particles comparable to the current experimental limit, 
and will be put to stringent tests as the next generation 
of experiments improves the sensitivity by another fac- 
tor of 100. Even though none of the experimental efforts 
have yielded positive results, the EDM searches, to this 
end, have ruled out many theoretical models. With ever- 
more-refined experimental techniques, the EDM searches 
continue to be of fundamental significance in particle and 
nuclear physics. 

The conventional technique used to measure EDM is 
the separated oscillatory fields method by Ramsey Q 
based on nuclear magnetic resonance: the EDM interac- 
tions induce an additional frequency shift in the Larmor 
precession, when the particle under study is subject to an 
electric field applied parallel (or anti-parallel) to a weak 
magnetic field. In the attempt to improve the experi- 
mental limit on the eEDM, we have been pursuing an al- 
ternative approach using a solid-state technique [6j . The 
application of a strong electric field to a paramagnetic 
insulator aligns the EDMs of valence electrons bound in 
the solid, leading to a small, yet non-zero magnetization. 
Even though the energy shift predicted from the EDM 
coupling for individual electrons is much smaller than 
the thermal energy, the cumulative effect from the large 
number of electrons in a solid sample leads to a net spin 
alignment equivalent of a few million Bohr magnetons. 
This Stark-induced magnetization can be detected using 
sensitive magnctometry. The use of Gadolinium Gallium 
Garnet (GGG, GdaGasO^) for the eEDM search was 
first proposed by Q and realized by § . The GGG ma- 
terial is advantageous for its high number density of the 
Gd 3+ ion (~T0 22 cm~ 3 ), each containing seven unpaired 
electrons on the 4/ shell, that leads to a strong magnetic 
response in a bulk sample. In addition, the GGG pos- 
sesses superb dielectric strength of 10 MV/cm and a high 
electrical resistivity that allows it to withstand a strong 
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electric field (necessary to perform eEDM measurements) 
with sufficiently small leakage currents. 

Ideally, the experiment has to be performed in an envi- 
ronment free of magnetic field, because the spin coupling 
to the magnetic field would certainly dominate over the 
small EDM interaction. In practice, even with the most 
hermetic magnetic shielding, some residual field is in- 
evitable. Therefore, the experiment is carried out in AC 
mode, in which the change of magnetization upon the 
reversal of the electric field is measured. Unfortunately, 
during the field reversal, transient currents create a mag- 
netic field that also flips direction. The transient field can 
die down quickly, but the presence of hysteresis effects 
that lead to finite remnant magnetization (with long re- 
laxation time) in the sample would be detrimental to the 
successful realization of this technique at the proposed 
sensitivity level. Therefore, we use a paramagnetic gar- 
net as a conservative measure against possible systematic 
effects arising from the use of ferromagnetic materials, 
even though the magnetic susceptibility is much higher 
in iron garnets. 
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FIG. 1: Volume magnetic susceptibility of the GGG sample as 
a function of temperature, measured at a maximum applied 
magnetic field of 400 Oe. The red curve shows the Curie- Weiss 
fit. 



II. SAMPLE CHARACTERIZATION 

The GGG belongs to the garnet family Q with a 
general structure of {^43j-[i?2](C3)0i2, where A denotes 
triply- ionized metallic ions, M 3+ , on a {dodccahcdral} 
lattice, B and C are ions on a [octahedral] and 
(tetrahedral) lattice, respectively. Oxygen ions, 2 ~, 
form a cage around the ions and balance the charge. In 
general, the couplings between these sub-lattices are anti- 
ferromagnetic (AF). The A, B, and C can be substituted 
by many different metallic elements, with varying degrees 
of magnetism, leading to a wealth of magnetic properties 
that change with temperature. The rare-earth elements, 
which are of interest to the eEDM search, can occupy 
the A sites. Due to relativistic effects, the EDM of para- 
magnetic atoms (ions) is enhanced by Z 3 , where Z is 
the atomic number. Thus, the EDM of the Gd 3+ ion 
(Z = 64) dominates over that of the Ga 3+ and 2 ~. 
Non-magnetic Ga 3+ ions populate the B and C sites, 
leaving the magnetic property of GGG to be determined 
solely by the spin interactions of the Gd 3+ ions on the do- 
decahedral lattice. Despite the intrinsic AF coupling, the 
spin of the Gd 3+ remains disordered and follows a typi- 
cal paramagnetic behavior with which the magnetic sus- 
ceptibility increases with decreasing temperatures. Only 
at temperatures lower than a few hundred milli-Kelvin 
does the AF coupling prevail, and the system becomes 
geometrically frustrated, transforming into a spin glass 
state [f(|. However, the spins are never fully ordered 
and the spin degree of freedom remains unfrozen at low 
temperatures [ll], [l2j . 

We measure the magnetic susceptibility x of polycrys- 
tallinc GGG samples (synthesized in our lab [|| using 
solid-state reactions [13[) using a SQUID susceptome- 
ter system 22[ . The volume magnetic susceptibility was 



measured at several values of the applied magnetic field 
(400 Oe, 100 Oe, and 10 Oe). The result with the 
maximum applied field of 400 Oe is shown in Fig. [TJ 
from 295 K to 2 K. The fit to the Curie- Weiss relation, 
X = C/(T — Ocw), gives the Curie- Weiss temperature of 
-2.1 K, indicating that coupling of adjacent Gd 3+ ions 
is indeed AF . The strong AF coupling could lead to an 
order-disorder phase transition at low temperatures, and 
limit the size of %■ Fortunately, for GGG, this phase 
transition is highly suppressed and was never observed 
due to the geometric frustration of AF coupled spins on 
a Kagome lattice [l(|. To maintain a high sensitivity to 
the eEDM, the spins need to remain free to respond to 
the external fields, therefore, it is essential to learn more 
about the conditions of phase transitions to ensure that 
the experiment is operated in the paramagnetic phase. 

Finally, to assess x °f the material, a correction for 
the demagnetization effect on the measured susceptibility 
Xmea is applied: 



X 



1 - N X r, 



(1) 



Here TV is the demagnetizing factor, which arises from an 
additional demagnetization field created by the magnetic 
surface charge density om = n-M. This leads to a partial 
cancellation of the applied field inside the sample, analo- 
gous to the electrostatic depolarization field. Our cube- 
shaped sample with dimensions of 0.3 cmx0.3 cmx0.3 cm 
is estimated to have N = 0.264, using finite-element- 
analysis calculations. 

We also measured the resistivity of the synthe- 
sized polycrystalline GGG sample using an electrome- 
ter (Kcithlcy 6517B), and attained a volume resistivity 
of (5.32 ± 0.04) x 10 15 f2-cm and a surface resistivity of 
(2.95 ± 0.02) x 10 15 Sl/cm 2 at 300 K. 
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gradiometer. It is placed between the ground electrodes. The 
inner diameter matches the diameter of the GGG sample. For 

FIG. 2: Schematic of the experimental setup for the GGG- best CMMR, the inner area irRj is matched to the outer area 

based solid-state eEDM experiment nR^ — -kR\. 



III. EXPERIMENTAL DETAILS 

With the samples characterized, we built the EDM ex- 
perimental cell from two disk-shaped GGG samples with 
diameter of 3.3 cm, height of 0.76 cm, and density of 
6.66 g/cm 3 , sandwiched between two planar electrodes 
and two isolated ground plates (Fig. [5]). The electrodes 
are connected to high voltage (HV) sources of oppo- 
site polarities in such a way that the electric fields in 
both GGG samples are aligned along the same direc- 
tion. In the presence of a strong electric field, the eEDM 
are aligned by the electric field, leading to a net spin- 
polarization, because the EDM vector is bound in the 
same (or opposite) direction of the spin vector as a re- 
sult of the Wigner-Eckart theorem. This Stark-induced 
spin ordering generates a bulk magnetization which pro- 
duces a magnetic field surrounding the paramagnetic 
GGG sample. The magnetization can be detectable us- 
ing a DC superconducting quantum interference device 
(SQUID), serving as a flux-to- voltage transducer. 

Due to the small energy of eEDM interactions, the in- 
duced magnetic field is minute. However, with EDMs 
near the current experimental limit, the accumulated 
magnetic signal from the large number of electrons in- 
side the GGG sample could lead to a sizable signal above 
the background. To be more explicit, with an internal 
field Ei n t of 10 kV/cm, the EDM interaction energy, 
d e Ei nt ~ 10 -23 eV, using the current upper limit for the 
eEDM, d e . In contrast, the thermal energy is as large 
as k B T = 8.6 x 10" 7 cV at 10 mK. where ks is the 
Boltzmann constant. The degree of EDM- induced spin 
alignment is washed out by the thermal fluctuation to the 
level of d e E mt /k B T ~ 10~ 16 for each Gd 3+ ion. On the 
other hand, the large number of Gd 3+ ions (each with a 
magnetic moment /z a ) in the solid results in a bulk mag- 
netization of 10 22 cm~ 3 xl0 _16 ^ a ~ 10 6 ^b, leading to a 
net magnetic field of 10~ 17 T. 

A flux pickup coil, in the form of a planar gradiome- 
ter, is sandwiched in between the ground planes. It in- 



tegrates the EDM-induced magnetic flux over the area 
of the coil. More importantly, the use of a gradiome- 
ter eliminates the common-mode magnetic signal from 
the residual magnetic field remaining inside the magnetic 
shield. It can significantly reduce the magnetic pickup 
due to the vibrational motions of the coil in a residual 
field. The gradiometer (Fig. [3]), with a two-turn inner 
coil wound clockwise and a single-turn outer coil wound 
counter-clockwise, is optimized to have the proper induc- 
tance matching to the input coil of the SQUID sensor. 
The diameter of the inner coil matches the diameter of 
the sample. The common-mode rejection ratio (CMMR) 
of a typical hand-wound coil was measured to be ~ 200, 
corresponding to a 0.5% area mismatch. The magnetic 
flux pickup is slightly enhanced by partially enclosing 
the returning flux, leading to an enhancement factor of 
1.1 (calculated using finite-element analysis) compared 
to that using a simple one-turn coil. In other words, the 
effective area A for flux pickup in this gradiometer is a 
factor of 1.1 higher than the actual cross-sectional area 
of the sample. This enhancement factor can be increased 
up to 1.8 by reducing the radial dimension of the super- 
conducting lead shield, and thus compressing the return 
flux lines laterally to increase the flux pickup. 

The EDM-induced magnetic flux in the sample en- 
closed by the pickup coil, $ e , can be estimated using [7j 



Ma 

Here a is the paramagnetic EDM enhancement factor for 
Gd 3+ ion in the GGG structure, which includes the effect 
of the dielectric reduction of the external field. E ext is 
the strength of the externally applied electric field, A 
is the effective area of the pickup coil aforementioned, 
/i Q is the magnetic moment of the paramagnetic Gd 3+ 
ion (p a = g^J J(J + l)/xg = 7.94//^, with g = 2 and 
J = 7/2), and / is the flux suppression factor due to the 
demagnetizing effects depending on the geometry of the 
GGG sample. 
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The effective EDM enhancement factor a in Eq. [2] was 
calculated by Sushkov's group [T^4l7j . In the perturba- 
tion calculation, the EDM-induced energy shift Ae per 
Gd 3+ ion arises from three independent effects, includ- 
ing (a) the EDM enhancement in the Gd 3+ ion, (b) the 
electron-electron coulomb interaction, and (c) the lattice 
deformation and the positional shift of the Gd 3+ ion with 
respect to the surrounding O 2- ions in the GdOs cluster. 
The resulting Ae is 35. 6d e Ei nt where Ei nt is the internal 
electric field inside the sample. Because of the dielectric 
reduction of the internal electric field of Ei nt = E ext IK ', 
where the dielectric constant K w 12 in GGG [181 fL9|. 
the overall energy shift per Gd 3+ ion is 

Ae = 35.64 f^f ) = 2.97d e E ext . (3) 

Note that the authors used a value of K = 30, quot- 
ing from an online table [23j. The references we found 
all suggest a smaller K ~ 12, and we also indepen- 
dently confirmed the smaller dielectric constant of the 
GGG with capacitance measurements. Hence, the effec- 
tive EDM enhancement factor a is 2.97, which is a factor 
of 3.5 larger than the original estimate, in which the en- 
ergy shift Ae is ~d a E l = - (-2.24)1 (K + 2)(E ext /K) 
= 0.86d e E ext , where d a = —2.2d e [l5j and Ei is a local 
field acting on Gd 3+ , estimated using a simple Lorentz 
relation Ei = \{K + 2)E lnt . 

The flux suppression factor / in Eq. [2] describes the de- 
gree to which the actual flux measurable from the EDM- 
induced magnetization is reduced due to the effect of 
geometry. The suppression factor / of the disk-shaped 
GGG sample with a 3.3 cm diameter and a 0.76 cm height 
is calculated using finite-element analysis. The solution 
shows a non-uniform magnetic field given a uniform mag- 
netization inside the disk-shaped sample. To obtain $ e , 
we integrate the solution field over the area of the pickup 
coil. The resulting magnetic flux is suppressed by a fac- 
tor of 0.369, which is a result of the finite dimensions of 
our sample geometry, and by another factor of 0.425 due 
to the placement of the pickup coil 0.33 cm away from 
the immediate surface of the sample. Note that this flux 
reduction is of a different origin from the demagnetiza- 
tion effect discussed before in the context of the magnetic 
susceptibility measurement. The total suppression factor 
/ is estimated to be 0.157, leading to a loss in sensitivity 
to the eEDM not considered in the original proposal [jj- 
Note that this suppression factor can be improved by 
moving the pickup coil closer to the sample. 

As shown in Fig. [21 the pickup coil (with an inductance 
L p of 618 nH) connects to the built-in input coil on the 
SQUID sensor chip (Supcracon CE2bluc, with an input 
coil inductance Li of 420 nH and a mutual inductance 
M between input coil to SQUID of 8.1 nH). The increase 
of magnetic flux 4> e picked up from the sample induces 
a current I = $ e /(L p + Li). The current flows into the 
input coil, and produces a flux $ sg that couples into the 
SQUID loop inductively, and is read out as a voltage 



signal. The relationship between $ s? and $ e is given by 

<P sq = MI = M $ e = /3$ e (4) 
Lp + Li 

where (3 is the coupling efficiency which quantifies how 
much the flux is diminished when <f> e is delivered to the 
SQUID sensor. The coupling efficiency is calculated to 
be 0.0078 in our setup. To enhance the measurable flux, 
we would need a strong electric field, a large sample size, 
and an optimized magnetic flux pickup coil (see Eq. [2]) . 

To reduce the Johnson noise, the electrodes are made 
of machinable ceramics (MACOR) coated with graphite 
to provide large but finite resistivity. The large electrical 
resistivity helps to reduce eddy currents and the mag- 
netic noise produced by random motion of conducting 
electrons. During the eEDM measurement, HV of oppo- 
site polarities are applied to the two HV electrodes so 
that the electric fields in both samples are in the same 
direction. The leakage current on each ground plate is 
monitored by a dedicated low-noise current preamplifier 
(Stanford Research Systems SR570). The assembly of 
the GGG samples and electrodes is shielded from exter- 
nal magnetic fields with two layers of superconducting 
lead foils and additional three layers of mu-metal (Met- 
glas alloy ribbon) wound on square forms with symmetry 
axes along the x, y, and z directions, respectively. The 
experiment is mounted inside a helium cryostat to al- 
low full immersion in a bath of liquid helium at 4.2 K at 
atmospheric pressure. Finally, a cylinder of Co-Netic fer- 
romagnetic shielding at room temperature encompasses 
the whole cryostat to provide the initial reduction of the 
ambient fields. The noise spectrum of the SQUID sensor 
instrumented inside the working experiment is shown in 
Fig. 2J Note that the experimental setup described here 
is the prototype design for the proof-of-principle mea- 
surement. Future work is planned to extend this work to 
sub-Kelvin temperatures, in order to attain better sensi- 
tivity through the enhanced magnetic susceptibility. 

IV. ESSENTIAL IMPROVEMENTS 

For stable operation of the SQUID sensor at the base 
noise level, any uncontrolled sources of electromagnetic 
interference (EMI) are undesirable. Considerable effort 
went into studying and eliminating electrical EMI noise 
in the lab, and eliminating ground loops. As shown in 
Fig. H] (dashed curve), the baseline of the SQUID sen- 
sor is 6/i<f>o/ \fB~z at the frequency of operation (close 
to the manufacturer's specification), after an improve- 
ment of over a factor of 50 in the magnetic shielding 
factor. The solid curve shows the spectrum with an ap- 
plication of a 2 Hz external uniform magnetic field of 
10 gauss generated by a pair of Hclmholtz coils placed 
outside the cryostat. By comparing the amplitude of the 
residual 2 Hz peak to the applied field strength, we esti- 
mate the overall magnetic shielding factor of the system 
to be 5 x 10 11 . The quality of the magnetic shield and 
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FIG. 4: The magnetic noise spectra of the DC SQUID detec- 
tor operated at 4.2 K. The nominal noise spectrum (dashed 
curve), and the spectrum with an external magnetic field of 
10 gauss at 2 Hz, applied outside the magnetic shield. 



FIG. 5: Time traces of averaged signals on one HV monitor, 
one leakage current monitor, and the SQUID readout (over 
three hours of sampling) presented as two cycles of field re- 
versal. Data were taken using the 24-bit DAQ system, with 
a sampling rate of 715 Hz. 



the ground-loop isolation also ensures that no observable 
vibrational peaks or power line harmonics show up on 
the SQUID noise spectrum. 

The expected EDM-induced magnetic signal to be 
measured by the SQUID sensor is small, thus any pos- 
sible contamination from other voltage monitoring chan- 
nels (such as the high voltage channel, which has a very 
large signal by comparison) through capacitive coupling 
is intolerable. To address this problem, a custom data ac- 
quisition (DAQ) system was designed, with the capability 
of eight-channel simultaneous sampling [2(|. The DAQ 
system has eight dedicated analog-to-digital converters 
(ADC) with 24-bit resolution for each analog input chan- 
nel, individually shielded in its own isolated heavy-duty 
radio frequency (RF) shielding enclosure. Fiber optic 
links are implemented between each satellite ADC board 
and the master board for measurement control and data 
retrieval. Using this DAQ system, we are able to mini- 
mize cross-talk between channels to better than -191 dB. 
The system is carefully designed to reduce electromag- 
netic interference, and eliminate the possibility of un- 
wanted currents flowing in the ground loops that result 
in increased noise. The intrinsic root-mean-square noise 
of the DAQ system is measured to be 1.38 /LtV. The DAQ 
system allows us to collect a large amount of data for av- 
eraging, without introducing additional sources of non- 
Gaussian noise at the level of the desired voltage sensi- 
tivity. The Gaussian-distributed random noise from the 
SQUID detector can thus can be reduced with higher 
statistics. 

This low-noise, high-resolution DAQ system enables us 
to study the correlations between the measured leakage 
current and applied high voltage, and uncover a system- 
atic effect originating from the voltage drift of the HV 
polarity switch system. After each polarity switch, the 
HV settling time could be quite long. With a small volt- 
age drift on top of the nominal 5 kV charging voltage, the 



unquenched charging and discharging of the electrodes 
leads to non-zero electric currents CdV/dt flowing in and 
out of the electrodes. This current is detected by the 
leakage current monitor. It also generates a magnetic 
field in phase with the polarity of the high voltage and 
leads to a signal (on the SQUID channel) that mimics 
the EDM-induced magnetization, given that the normal 
of the pick-up coil is not perfectly aligned with the field 
lines. To eliminate this source of background, we im- 
proved our HV polarity switch system to have distortion 
<10 ppm/sec of voltage drift. In order to handle the po- 
larity switching between ±5 kV, we use vacuum-tube tri- 
odes (6BK4C) connected in series to the positive and neg- 
ative DC HV supplies (Stanford Research PS350). The 
gate voltage for each vacuum tube is controlled through 
a opto-isolator, driven by an arbitrary waveform genera- 
tor; the input square waveform can be amplified by the 
HV polarity switch system by a factor of 1000. A re- 
duction of the voltage drift by a factor of 600 over the 
previous supply was achieved by improving the feedback 
circuit and reducing the transient ramp at field reversals. 
Currently, the voltage drift on the HV output is limited 
by the drift of the low-voltage drive from the arbitrary 
function generator which has only 14-bit resolution. Our 
next-generation HV system uses a precision 22-bit DAC 
drive (integrated to the DAQ system) and should im- 
prove the dV/dt by another factor of 30, with increased 
HV output from 10 kV PJ3 to 40 kV pp . 

In addition, the early problems we encountered with 
SQUID instabilities, that resulted in frequent flux jumps 
at applied electric fields higher than 3 kV/cm, were re- 
solved by replacing HV cables and adding low-pass filter 
resistors, which needed to be cryogcnically-compatiblc 
and HV-rated. This improvement allows the application 
of the full range of 10 kV pp across the sample. 
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Fit model: y=Ae (x/t) +Bx+d 



FIG. 6: Distribution of the eEDM observable (processed with 
the drift-correction algorithm). Data points were collected 
over one non-interrupted run (3 hours). 
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FIG. 7: SQUID readout (averaged over one non-interrupted 
run) folded modulo one cycle of field reversal. This data set 
illustrates the worst case of DC voltage drift. 



V. ANALYSIS 

Typical averaged time-traces of three hours of data 
on the monitored channels are plotted in Fig. [5j In an 
eEDM measurement sequence, we applied a voltage of 
alternating polarities up to ^10 kV pp in a square wave- 
form with controlled ramp speed on the electrodes. The 
polarity switching cycle is repeated at a rate of 1.43 Hz. 
This drive frequency is chosen to be low enough to re- 
duce the transient currents, but high enough to avoid the 
low frequency 1// corner of the SQUID noise spectrum 
(Fig.[4j). Using the 24-bit DAQ system, we monitored the 
high voltage on the two electrodes (through 1000:1 volt- 
age dividers), currents flowing in the two isolated ground 
plates, and the analog voltage output from the SQUID 
readout electronics that could be converted into the mag- 
netic flux $ sg through a predetermined transfer function. 

The time-trace of the current monitor contains both 
the charging and discharging current, CdV/dt, and the 
leakage current, V/R, flowing through the bulk sample 
or the surface. As shown in Fig. [5J the measured cur- 
rents through the ground plates are dominated by the 
charging/discharging transient currents during the HV 
polarity switching. The SQUID sensor measures mag- 
netic fields, generated by the EDM-induced sample mag- 
netization as well as by the electric currents flowing in 
and out of the electrodes. During field reversals, the 
SQUID is measuring the large magnetic field associated 
with transient currents, which have to die down quickly 
so that the SQUID can measure the EDM-induced mag- 
netic flux once the field settles to the maximum ampli- 
tude of the applied HV. The difference of the magnetic 
fluxes at +HV and -HV within one cycle is proportional 
to the eEDM, and is defined as our eEDM observable. 
The dV/dt needs to be controlled during the time win- 
dow of eEDM measurement (last 30% of the half cycle) , 
so as not to dominate over the magnetic flux generated 
by the induced magnetization. 



Despite all the improvements discussed above, there 
still exist residual voltage drifts. The worst case of such 
drift is illustrated in Fig. The DC drift could come 
from many sources, including the SQUID electronics, the 
slow reduction in the level of the liquid helium, and the 
pressure drift inside the cryostat, to name a few. Unlike 
the drift in the HV source, which changes sign in phase 
with the HV polarity, this DC drift does not have the 
same correlation with the HV cycle. To remove this DC 
drift from the data, we use two independent algorithms 
(drift-correction and fitting) to analyze the eEDM data. 
The get the eEDM observable, the drift-correction algo- 
rithm takes the algebraic sum of the SQUID readout in 
two adjacent cycles of field reversal, and applies a [+1 -3 
+3 -1] weighting on the averaged data for each half cy- 
cle. With this, the effect of the DC drift can expand into 
a polynomial function of the time, and canceled up to 
second-order. The transient regions at the field reversal 
are excluded in the data average. Furthermore, to ensure 
that the transient current has sufficiently decayed, the 
data window contains only the last 30% of the time trace 
(Fig. [5]). Fig. [5] shows the histogram of eEDM observable, 
collected over a typical run with three hours of data. The 
distribution can be fitted by a Gaussian distribution, giv- 
ing <E> S? = (-1.89 ± 9.25) x 10" 8 $o- This corresponds to 
an eEDM measurement of (-0.24 ± 1.17) x 10~ 23 e-cm, 
after taking into account of all the suppression factors 
discussed in the previous section. 

The fitting algorithm attempts to fit the entire time 
trace of the SQUID readout modulo one cycle of field 
reversal. Data from repeated cycles over three hours are 
averaged to form the overall time-trace (Fig. [7]), and then 
fitted using the following voltage model: 

V + {t)=Ae- t ' T +Bt + d + , 

V- (t) = -A'e-V r ' + {Bt + B^) + d~, 

where V + (t) and V~(t) are the SQUID readout during 
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5x10" 
4x10" 
3x10" 
2x10" 
1x10" 



TABLE I: Systematic effects 



'-1x10" 
-2x10" 
-3x10" 
-4x10" 
-5x10" 



Individual d e extracted by the fitting-algorithm 
Final d 




(-0.77+/-7.98)x10 25 e-cm 



5 10 15 20 25 30 35 40 
Run number 



Source 




Leakage current 
Displacement current 
Remnant magnetization 
Channel cross-talks 
Vibrational peak © 2.13 Hz 


(3.18 ± 0.07) x 10- a $ ° 
1.65 x 10- 10 *o b 

< 6.02 x 10~ 9 <I>o c 
(0.52 ± 1.51) x 10~ 10 $ d 

< 6.02 x 10" 9 <E>o e 


Current Flux Measurement 


(0.48 ±6.02) x 10~ u $ o ; 



"Average over all runs. 
b C^- < 0.2 pA at 9.4 kV pp . 

^Direct MPMS measurements are limited by the remnant fields. 
d Channel isolation> 191 dB. 

e No observable peak on the SQUID noise spectra. 
^Averaged over all runs. 



FIG. 8: Results of eEDM measurements, processed using the 
fitting algorithm. Each data point is the averaged result of 
one run. The final data point is the final average results over 
40 experimental runs. 

the half-cycle with either a positive or a negative HV 
polarity. The first term in the model characterizes the 
decay of the transient current with a time constant r 
and t'. These two time constants could be different be- 
cause of the asymmetry of the circuit handling the posi- 
tive and negative voltage, the difference between the two 
HV channels, and/or the two HV electrodes. The sec- 
ond term describes the DC voltage drift. It is Bt for the 
first half cycle, and Bt + B(T/2) for the second half a 
cycle, which starts (T/2) later, where T is the period of 
the field reversal cycle. The last constant term repre- 
sents the EDM- induced magnetization and an DC offset. 
Note that the sign of the EDM-induced magnetization 
changes as the electric field is reversed, while the DC off- 
set remains constant. The eEDM observable is derived 
simply by taking the difference of the fitted parameters 
d + and d~ , 

AV = d + -d~. (6) 

The results of the fitting algorithm give $ sq = (3.07 ± 
6.34) x 10™ 8 $o, corresponding to a d e of (0.39 ± 0.81) x 
10 -23 e-cm. The fitting algorithm arrives at a better 
statistical sensitivity than the drift-correction algorithm 
because the fitting algorithm uses about 75% of the col- 
lected data points, as opposed to 30% used in the drift- 
correction algorithm. We collected data over two weeks, 
with a total integration time of five days. 

Analysis using the two algorithms shows similar dis- 
tributions of the extracted eEDM, and the results us- 
ing the fitting-algorithm are shown in Fig. [3] Each data 
point corresponds to an non-interrupted run lasting for 
three to four hours. Proper error-weighting is applied to 
sum the results of each run to arrive at a final eEDM of 
(0.41±1.38) x 10~ 24 e-cm and (-0.77±7.98) x 10~ 25 e-cm, 
for the drift-correction and fitting algorithms, respec- 
tively. These results can be compared to the previ- 



ous experimental limit using a complimentary solid-state 
method in the gadolinium iron garnet system pll ] . 



VI. DISCUSSION 

A comprehensive list of systematic effects is shown in 
Table Q] Since the measured physical observable is the 
magnetic flux, we compare the spurious flux generated 
by each known systematic effect. The dominant effect is 
the leakage current through the sample which produces 
a magnetic field in phase with the polarity of the HV. 
To first order, the magnetic field generated by the leak- 
age current is perpendicular to the EDM-induced mag- 
netization, some fraction of the field, however, can be 
measured by the SQUID sensor due to the slight tilt 
of the pickup coil. Surface currents forming a helical 
path would generate additional magnetic flux. Studies of 
the correlation between the displacement currents (dur- 
ing field transients), the time-derivatives of the applied 
HV, and the SQUID signal allow us to separate the con- 
tributions from the displacement current from that of 
the leakage current (due to the finite resistivity of the 
sample). These studies show that about 1.4% of the ra- 
dial field generated by the the leakage current could leak 
into the pickup coil and contribute to a spurious signal. 
For example, a leakage current of (6 ± 2) pA (averaged 
over three hours) measured at the maximum applied volt- 
age of 10 kVpp generates a spurious magnetic flux $ s? of 
(3.79 ± 1.26) x 10 _9 $o- Despite the finite leakage cur- 
rent, the measured eEDM is shown to be independent of 
the strength of the applied electric fields within the error 
bars, as shown in Fig. [9] Most measurements were made 
at 9.4 kVpp. This suggests that the experiment is free of 
systematic effects linear in the HV. 

We also attempted to characterize the magnetic hys- 
teresis of the GGG sample using the SQUID-based sus- 
ceptometer system, however, due to the large variations 
of the remnant field upon each field ramping, these mea- 
surements cannot be used to determine the magnetic hys- 
teresis of the GGG sample to the sensitivity level required 
in our experiment. Since there is no measurable offset be- 
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FIG. 9: Measured eEDM as a function of the applied HV. 
The data point at 9.4 kV pp contains the highest statistics. 



tween the half cycles, we can place an upper limit on the 
remnant magnetization using our EDM results. Given 
the knowledge of the magnetic susceptibility of GGG, 
we need to control the remnant magnetization, which 
changes in phase with the applied field, to be smaller 
than 2.5 x lOVs/cm 3 . 

The dominant systematic effect is the leakage current, 
which creates an additional magnetic flux, and can lead 
to a spurious eEDM signal through Eq. [5J The mea- 
sured flux is corrected for each run by subtracting the 
additional magnetic flux created by the leakage currents, 
before extracting the eEDM value. The total systematic 
effect is estimated to be (4.80±0.12)x 10~ 25 e-cm with the 
proper error-weighting. This leads to the final reported 
eEDM value of (-5.57±7.98 stQ t±0.12 syst ) x 10" 25 e-cm 
with five days of data averaging, in this prototype exper- 
iment running at 4.2 K. 



temperatures using a dilution refrigerator. At such tem- 
peratures, the spin ordering in the GGG sample would 
increase according to a 1/T behavior of the magnetic 
susceptibility, so that the EDM sensitivity is expected to 
increase by about a factor of 100 at 40 mK, without any 
improvements to other parts of the experiment. However, 
the non-zero Curie- Weiss temperature and the spin-glass 
phase transition could potentially limit this gain. It will 
require careful measurements of the magnetic properties 
of the GGG material at sub-Kelvin temperatures to as- 
sess the impact of its intrinsic AF spin coupling. The 
use of the garnet allows us to perform further systematic 
studies by substituting the magnetic Gd 3+ ions by non- 
magnetic Y 3+ ions on the same lattice sites. Our prelim- 
inary tests have shown that the spin dilution changes the 
Curie- Weiss temperature, and it could also potentially 
help by pushing the spin-glass phase transition to lower 
temperatures. 

The expected improvement in eEDM sensitivity origi- 
nates from the enhanced magnetic susceptibility at sub- 
Kelvin temperatures. It requires same level of uncertain- 
ties in the flux measurements, and thus the requirement 
of control on all the systematic effects discussed above 
do not need to be further improved. In fact, the leading 
systematic effect from the leakage current can only be re- 
duced by cooling the experiment to lower temperatures, 
at which the resistivity of the sample would increase. In 
addition to low-temperature operation, we plan to scale 
up the size of the experiment by using 10 modules of 
the sample/electrode assembly tested in this prototype 
experiment, increase electric fields up to 40 kV pp , and 
employ a better SQUID sensor with lower noise and in- 
creased coupling efficiency. With these planned improve- 
ments, we expect to improve the experimental sensitiv- 
ity of the solid-state technique to and beyond that of the 
atomic beam experiments [4J. 



VII. CONCLUSIONS AND FUTURE WORK 

In this work, we have learned to solve many systematic 
problems, and demonstrated the feasibility of the solid- 
state method using the GGG paramagnetic insulator at 
4.2 K. We report our first background- free experimen- 
tal limit on the eEDM. Further improvement of the ex- 
perimental sensitivity using the GGG system is possible. 
It would require cooling the experiment to milli-Kelvin 
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